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Abstract 

A Chern insulator (quantum anomalous Hall insulator) phase is demonstrated to exist in a typical 
semi-Dirac system, the Ti02/V02 heterostructure. By combining first-principles calculations with 
Wannier-based tight-binding model, we calculate the Berry curvature distribution, finding a Chern 
number of —2 for the valence bands, and demonstrate the existence of gapless chiral edge states, 
ensuring quantization of the Hall conductivity to 2e^ jh. A new semi-Dirac model, where each semi- 
Dirac cone is formed by merging three conventional Dirac points, is proposed to reveal how the 
nontrivial topology with finite Chern number is compatible with a semi-Dirac electronic spectrum. 

PACS numbers: 73.20.-r,73.21.-b,73.43.-f 
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A Chern insulator is a two-dimensional (2D) magnetic insulator with a quantized anoma¬ 
lous Hall conductivity Ce^/h, where C is an integer topological index known as the Chern 
number These systems, also known as quantum anomalous Hall (QAH) insulators, have 
attracted a great deal of interest, in part because of their gapless chiral edge states which 
enable dissipationless transport, with potential applications in electronic devices j^. So 


arjSeveral systems have been proposed, notably magnetically doped topolqgica’ 


lattices 


12 


, 1^ leading to recent experimental conhrmation j^, but also honeycomb jbl-ll l| or square 


r 


insulators 


14| formed by transition-metal and heavy-metal ions. The essential ingredients 


are the spontaneous breaking of time-reversal (TR) symmetry, as by the formation of a 
ferromagnetic state, and the presence of spin-orbit coupling (SOC), which generates the net 
Berry curvature needed for a nonzero Chern number. Here we focus on one interesting class 
of proposals involving 2D systems that would be half-semimetals in the absence of SOC, 
with the Fermi energy pinned at one or more Dirac points in one spin channel while the 
other channel is gapped. If the application of SOC gaps the Dirac points to produce conical 
avoided crossings instead, a nonzero Chern number can result, as proposed for example for 
a triphenyl-manganese (Mn 2 Ci 8 Hi 2 ) system {q]. 

Recently, calculations on Ti02/V02 multilayer structures led to a proposal jlA. 16| for a 
new kind of half-semimetal in which, in the absence of SOC, the bands in the ungapped spin 
channel have a semi-Dirac dispersion, i.e., a quadratic rather than a linear dispersion in one 
direction. This raises the interesting question whether such a system can also provides a 
route to a Chern-insulating state when SOC is included. 


In this work we answer this question. First, we distinguish between two qualitatively 


differen 
Ref. 


types of semi-Dirac structures. The one we denote as “type-I” was proposed in 


17j . but we hnd that it cannot lead to a QAH state. Instead, we hnd that the Ti02/V02 
multilayer structure [^, ^ is described by a new “type-H” semi-Dirac cone structure, which 
does lead to a Chern-insulating state when SOC is turned on. Moreover, we clarify that the 
type-H structure is not protected by symmetry, and (in the absence of SOC) will generically 
transform into one or three Dirac nodes in the absence of hne tuning. With SOC, we 
predict that the Ti02/V02 heterostructure is a QAH insulator with a Chern number of —2, 
demonstrating a new route to the formation of a Chern-insulating state in 2D. 
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Dirac and semi-Dirac cones. A general effective two-band Hamiltonian in 2D can be 
given as 

H(k) = h(k) • a, (1) 

where a = o'y,az) are Pauli matrices. In 2D k-space a form like 

h{k) = {vFk:,;,VFky,0) ( 2 ) 


describes a massless Dirac cone structure, appropriate to the case that SOC is absent, 
where wp is the Fermi velocity. SOC is taken into account by adding a mass term such that 
h(k) = {vFkx,VFky,mz), opening an energy gap of 21771^1 at the Dirac point. The Chern 
number is determined by integrating the Berry curvature D(k) over the 2D Brillouin zone 
(BZ), but for weak SOC the dominant contributions will come from Berry fluxes $ = ±7r 
concentrated near the avoided crossings. This follows from the well-known form of the Berry 
curvature for the two-band model of Eq. ([1]), which yields 


$ = / dk D(k) = / dk 


(9h 9h, 

X ^). 


2|hP dk^ dk. 


(3) 


y 


Hence each gapped (massive) Dirac cone makes a contribution of ±| to the total Chern 
number C, which is therefore determined by summing over all the massive Dirac cones at 
the Fermi level. 

On the other hand, a semi-Dirac spectrum is a peculiar energy dispersion in which 
quasiparticles behave as massless along one principal axis but as massive fermions along 
the perpendicular direction. Recently this novel spectrum was observed in a multilayer 
(Ti02)m/(V(Wn nanostructure (m > 5;n = 3 or 4) by Pardo et al. using hrst-principles 
calculations [l5|, . They proposed that the semi-Dirac electronic spectrum could be de¬ 

scribed by Eq. ([I]) with 17 


19| 


h(k)= (4) 

or a similar expression. We refer to this model as a “type-I” semi-Dirac model to distinguish 
it from the “type-H” model we propose below. Obviously the dispersion is massless (linear) 
along the ky direction but massive (quadratic) along k^, satisfying the dehnition of a semi- 
Dirac cone. Since Eq. (|1]) is an even function of it is evident that the Berry curvature 
D(k) is an odd function of kx according to Eq. ([3]). Hence the Chern number, given by 
the integral of D(k) over the vicinity of the semi-Dirac point, is zero after adding a mass 
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term H' = rUzCz to open an energy gap. This is consistent with the observation that the 
semi-Dirac point under Eq. (|1]) carries a zero Berry flux [^, in contrast with conventional 
Dirac points which provide a Berry flux of ±7r. 

We therefore conclude that type-I semi-Dirac cones have a trivial topology with Chern 
number C = 0 and thus do not lead to a Chern-insulating state. Nevertheless, by using 
first-principles calculations and a Wannier-based tight-binding analysis, we And that SOC 
does turn the (Ti02)m/(V02)n into a Chern insulator. We first present the results and then 
explain them. 

Methods. The first-principles electronic structure calculations are performed within the 
framework of density functional theory (DFT) 20|] as implemented in the QUANTUM- 


ESPRESSO package |2l| with the plane-wave pseudopotential method, and in the VASP 


package 


22l | with the projector augmented-wave method. We adopted the Perdew-Burke- 


Ernzerhof generalized gradient approximation (GGA) exchange-correlation functional 23 |. 
The kinetic energy cutoff is fixed to be 500 eV and a P-centered 8x8x1 k-point mesh is used 
in all cases. We treat the open-shell 3d orbitals by adding an effective Hubbard U correction 
of 3.4 eV on V and Ti atoms within the GGA-I-U approach 2J, l25|. Structural relaxations 
are carried out without SOC, and then the electronic structure is computed twice, once 
with and once without SOC. Ferromagnetic ordering is found to be energetically favored 
in both cases, and the ordering is along the z axis when SOC is present. We use Wannier 
interpolation based on maximally localized Wannier functions (MLWFs) to calculate the 
Berry curvature and the anomalous Hall conductivity, which requires a very dense k-point 


grid in the BZ 


26, 3- 


Chern insulator behavior of the (Ti02)5/(V02)^ system. We focus on the (Ti02)m/(V02) 
system with n = 3 and m = 5; because Ti02 is strongly insulating, the latter is enough to 
effectively separate the VO 2 trilayers into isolated 2D systems, in which the energy bands 
close to the Fermi level are dominated by the 3d states of V. Since there are two kinds of 
VOe octahedra whose in-plane O-V-0 chains are perpendicular to each other, we align our 
local axes differently on the two V sites as shown in Fig. [U^a). Because of the distortion of 
the VOg octahedra away from cubic symmetry, the triply degenerate t 2 g orbtials of V ions 
split into two doubly degenerate d± orbtitals {yz and xz) and one d\\ orbital — y"^) p 


We projected the Bloch wavefunctions onto these local orbitals to get MLWFs, which serve 
as an ideal basis for further analyses. As shown in Fig. [H The MLWFs keep the shape and 
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FIG. 1. (Color online) Amplitude isosurfaces of atom-centered d-like maximally localized Wannier 
functions (red for positive values and blue for negative), (a) d_\_ {yz and xz) and (b) dy (x^ — y^) at 
different V sites in the (Ti02)5/(V02)3 superlattice. Local coordinate systems are shown in green; 
local z and y axes are along the in-plane O-V-0 chain and the global c axis respectively. 


symmetry of local atomic d orbitals. 

Based on these atom-centered d-like MLWFs, we calculate the Wannier-interpolated en¬ 
ergy bands close to the Fermi level 29(], shown in the top panel of Fig. [2](a). The Wannier- 
interpolated energy bands are in excellent agreement with the DFT results (not shown). 
The semi-Dirac character is clearly apparent, with two bands crossing linearly along the 
F-M line and quadratically along A-B (perpendicular to F-M). When SOC is considered in 
our calculation, an avoided crossing occurs and a small band gap of ~2.5 meV opens at the 
crossing point. To investigate the topological nature of this system, we calculate the Berry 
curvature D(k) of all states below the Fermi level using 301 

(t/^nk | V^rnk) (t/^mk | t/^nk) 


^(k) = - 2Im- 

n<EF m^n 




(5) 


where V'nk is the spinor Bloch wavefunction of band n with corresponding eigenenergy 
and V = {vx,Vy) is the velocity operator. In the bottom panel of Fig. EKa) the Berry 
curvature is plotted along the same /j-path. The large peaks between the F and M points 
arise where the conduction and valence bands are nearly degenerate and only weakly split 
by SOC, giving rise to small denominators in Eq. ([5]) and hence a large contribution to r2(k). 
The peaks show strong anisotropy, being much sharper along F-M than along A-B. The 
2D plot in Fig. [2](b) makes it clear that these avoided crossings give rise to four banana- 
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FIG. 2. (Color online) (a) Band structure and Berry curvature —II(k) along and perpendicular 
to the diagonal of the BZ. (b) Berry curvature —n(k) in the entire BZ. (c) Anomalous Hall 
conductivity axy plotted with respect to the position of the Fermi energy Ei?. The inset shows axy 
with E^;’ varying from —400 to 400 meV. (d) Energy and momentum dependence of the LDOS at 
the edge along the [11] direction. Note that the color maps are logarithmic; deep blue patches in 
(b) are numerical artifacts. 


shaped peaks of Berry curvature located around the semi-Dirac points. Note that f2(k) is 
not an odd function along the A-B direction, in contrast with the prediction of Eq. (|4j). 
Since the four peaks are related by fourfold rotations, there is no cancellation among them; 
each contributes a Berry flux of — tt, so that the total Chern number is —2. The system is 
therefore a Chern insulator. 


We also plot the intrinsic anomalous Hall conductivity [311 as a function of Fermi energy 
in Fig. |2][c). As expected from the nonzero Chern number, the anomalous Hall conductivity 
shows a quantized Hall plateau at axy = —Ce^/h when the Fermi level lies inside the bulk 
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band gap. The width of the Hall plateau is about 2.5 meV, corresponding to the nontrivial 
bulk band gap. The anomalous Hall conductivity decreases rapidly to zero when the Fermi 
level is outside the band gap, as expected since only states near the nontrivial bulk band 
gap, which is opened by the weak SOC, contribute strongly to the Berry curvature [see Fig. 
2(a,b)]. 

The existence of topologically protected chiral edge states is one of the most important 
consequences of the QAH state. To further reveal the nontrivial topological nature of the 
system, we calculate the edge states of a semi-inhnite (Ti 02 ) 5 /(V 02)3 system with its edge 
along the [11] direction. Because the existence of chiral edge states is completely determined 
by the bulk topology, here we build the tight-binding model with SOC based on first- 
principles MLWFs from the bulk, ignoring the effects of edge reconstruction. We apply an 


iterative method 


32| to obtain the edge Green’s function and the local density of states. 


which is directly related to the imaginary part of the Green’s function. In Fig. [2][d) we can 
clearly see that two edges states connect the valence and conduction bands. The appearance 
of two chiral edge states is consistent with our calculated Chern number C = —2, confirming 
the nontrivial topological nature of this semi-Dirac system. 

Type-II semi-Dirac model. The above results demonstrate that the (Ti02)5/(V02)3 su¬ 
perlattice is undoubtedly a Chern insulator, contrary to previous understandings that a 
semi-Dirac cone structure should not result in a Chern-insulating state. To resolve such a 
conflict, we propose a “type-H” semi-Dirac model with an effective 2x2 Hamiltonian in the 
absence of SOC having the form of Eq. ([T]) with 


/ 

h(k) = ( ^ otkxky, 0 


( 6 ) 


or other equivalent expressions. The resulting dispersion relation satisfies the semi-Dirac 
character: it is linear along one direction [k^ = 0) and quadratic perpendicular to this 
direction [ky = 0). Hence the origin point (0,0) is a semi-Dirac point. We then add a 
SOC-induced mass term H' = m^Cz to open a gap at the semi-Dirac point. 

To understand the nature of the Berry curvature expected from the type-I and type-H 
models of Eqs. dl]) and ([H]), we plot band dispersions and Berry curvatures derived from 
these models in Fig. [21 In particular, we investigate the evolution of the crossing points 
when hx aquires a small perturbation A in both models (assuming m > 0 without loss of 
generality). In the type-I model this leads to the the fusion of two Dirac points as A passes 
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FIG. 3. (Color online) Evolution of band structures of (a)-(c) type-I and (d)-(f) type-II semi-Dirac 
dispersions when hx aquires a small perturbation A. Berry curvature —D of (g) type-I and (h) 
type-II semi-Dirac models with A = 0 when a SOC-induced mass term H’ = mzOz is added to 
open a small gap at the semi-Dirac point. 

through zero. As shown in Figs. 3(a-c), when A < 0 there are two Dirac cones located 
at k = (±\/—2mA, 0); these merge to form a semi-Dirac cone at A = 0, and then a gap 
of magnitude 2A appears for A > 0. Because the two Dirac points provide Berry fluxes 
of opposite sign, the total Berry flux vanishes, as shown for A = 0 in Fig. [3](g). Thus, the 
type-I model is topologically trivial, and cannot be used to describe the semi-Dirac structure 
observed for the (Ti 02 ) 5 /(V 02)3 superlattice system. 

For the type-II model, on the other hand, when A < 0 the energy gap vanishes at three 
Dirac points located at k = (0, A/ui?) and (±\/—2mA, 0), as shown in Fig. [3](d). As A 
increases, these approach each other and merge at A = 0 into a single semi-Dirac point at 
k = (0, 0) [Fig. [3Ke)]. For A > 0, only the touching at k = (0, A/vp) remains; the other two 
points disappear (their coordinates formally become imaginary) [Fig. [3][f)]. Therefore, the 
type-II semi-Dirac model of Eq. ([6]) can be viewed as the consequence of the merging of three 
Dirac points to produce a single Dirac point. The distribution of the Berry flux produced by 





the type-II model at A = 0 is shown in Fig. [3]^h). Signihcantly, the Berry cnrvatnre near the 
semi-Dirac point has the banana shape shown in Fig. [3l^h), very similar to the DFT results 
in Fig. [2]^b). 

For one sign of ruz, the total Berry flux contributed by the type-II model is -I-tt indepen¬ 
dent of A; the two Dirac points at (±\/—2mA, 0) contribute +ti when A < 0, while the one 
at k = (0, A/uir) switches from —tt to as A crosses through zero. For the other sign of 
m^, all these are reversed. Since there are four symmetry-related semi-Dirac points in the 
entire BZ, the total Chern number is ±2. The results of the type-II semi-Dirac model thus 
agree very well with the DFT calculations of (Ti02)5/(V02)3 above, strongly suggesting 
that the type-II model correctly describes this system. 

To further verify this conclusion, we conduct an effective k ■ p analysis on the first- 
principles data without SOC around the band crossing point at (fcc, kc) and introduce a 
coordinate transformation q = (gi, qA with qi = {kx + ky)/2 — kc and q 2 = {kx — ky)/2. By 


employing the downfolding method 


33, 3 


we obtain the 2x2 effective Hamiltonian 35 1 


-ffk-p = e(q)l2x2 + h(q) ■ (T, 


(7) 


where 12x2 is the unit matrix. Keeping only terms up to quadratic order in q we find 
h(q) = {Aqi + Bq 2 , Cq 2 + Dqiq 2 , 0). The coefficient C is small but non-zero, implying that 
the dispersion along q 2 is almost quadratic, but contains a small linear component. We 
have checked that the presence of a nonzero coefficient C is not prevented by any symmetry. 
Setting h = 0 we hnd three zero-gap Dirac points: a real one at gi = g 2 = 0 and two virtual 
ones at (0.20, ±0.56i) (in units of a~^, where a is the in-plane lattice constant). These results 
map onto the A-modihed type-II model as discussed following Eq. ([6]), with qi = ky — A/vp 
and q 2 = kx, so we can identify A = —vp, B = l/2m, C = aA/vp, and D = a. Putting in 
numbers, we hnd a small positive A, placing the system into the regime of Fig. [3](f), but close 
to Fig. [3](e). If we artihcally decrease C and let it pass through zero, the real and virtual 
Dirac points merge at q = (0, 0) and reemerge as three real ones, exactly as predicted in 
Fig. [3](d) for the type-II semi-Dirac model. 

We therefore conclude that the (Ti02)5/(V02)3 nanostructure is an anisotropic linear 
Dirac system that is so close to a type-II semi-Dirac behavior that it can hardly be distin¬ 
guished from it. When the system is gapped by SOC, the calculated Berry hux is $ ~ —tt 
and the Berry curvatures are banana-shaped, independent of whether the linear term of q 2 is 
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exactly zero or not. 3^ Hence the Chern insulator phase is quite robust, even if the spectrum 
diverge from the exact semi-Dirac dispersion due to external strain or other perturbations. 


Conclusion. We have found that the (Ti02)5/(V02)3 nanostructure, which exhibits a 
peculiar semi-Dirac-like electronic structure in the absence of SOC, is a QAH (Chern) in¬ 
sulator with C = —2 when SOC is included. The calculated Berry curvature, anomalous 
Hall conductivity and gapless edge states unambiguously identify the nontrivial topological 
nature of this system. We propose a new semi-Dirac model, which describes the semi-Dirac 
point as the merging of three conventional Dirac points, to explain how a Chern-insulating 
state is compatible with a semi-Dirac dispersion. Our results clarify when and how the 
QAH state can emerge from a semi-Dirac structure, providing a new potential route to the 
realization of such states in practical materials systems. 
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Note added 
Cai et al. 


- Recently, a QAH phase was also found in Ti02/Cr02 heterostructures by 
. However, it exhibits nearly isotropic Dirac cones, which is quite different 


from the semi-Dirac cones of Ti02/V02 in our work. 


[1] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

[2] B. I. Halperin, Phys. Rev. B 25, 2185 (1982). 

[3] C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang, Phys. Rev. Lett. 101, 146802 (2008). 

[4] R. Yu, W. Zhang, H.-J. Zhang, S.-C. Zhang, X. Dai, and Z. Fang, Science 329, 61 (2010). 

[5] C.-Z. Chang and et al., Science 340, 167 (2013). 

[6] Z. Qiao, S. A. Yang, W. Feng, W.-K. Tse, J. Ding, Y. Yao, J. Wang, and Q. Niu, 
Phys. Rev. B 82, 161414 (2010). 

[7] D. Xiao, W. Zhu, Y. Ran, N. Nagaosa, and S. Okamoto, Nat. Commun. 2, 596 (2011). 

[8] H. Zhang, C. Lazo, S. Bliigel, S. Heinze, and Y. Mokrousov, 

Phys. Rev. Lett. 108, 056802 (2012). 


10 












[9] Z. F. Wang, Z. Liu, and F. Liu, Phys. Rev. Lett. 110, 196801 (2013); H. Hu, Z. Wang, and 
F. Liu, Nanoscale Res. Lett. 9, 690 (2014). 

[10] K. F. Garrity and D. Vanderbilt, Phys. Rev. Lett. 110, 116802 (2013), 

[11] G.-F. Zhang, Y. Li, and G. Wu, Phys. Rev. B 90, 075114 (2014). 

[12] H. Zhang, J. Wang, G. Xu, Y. Xu, and S.-C. Zhang, Phys. Rev. Lett. 112, 096804 (2014), 

[13] K. F. Garrity and D. Vanderbilt, Phys. Rev. B 90, 121103 (2014). 

[14] H. Zhang, H. Huang, K. Haule, and D. Vanderbilt, Phys. Rev. B 90, 165143 (2014), 

[15] V. Pardo and W. E. Pickett, Phys. Rev. Lett. 102, 166803 (2009), 

[16] V. Pardo and W. E. Pickett, Phys. Rev. B 81, 035111 (2010). 

[17] S. Banerjee, R. R. P. Singh, V. Pardo, and W. E. Pickett, 

Phys. Rev. Lett. 103, 016402 (2009), 

[18] G. Montambaux, E. Piechon, J.-N. Euchs, and M. O. Goerbig, 

Phys. Rev. B 80, 153412 (2009); G. Montambaux, F. Piechon, J.-N. Fuchs, and M. O. 
Goerbig, Eur. Phys. J. B 72, 509 (2009). 

[19] S. Banerjee and W. E. Pickett, Phys. Rev. B 86, 075124 (2012), 

[20] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964), 

[21] P. Giannozzi and et al., J. Phys.: Gondens. Matter 21, 395502 (2009), 
http://www.quantum-espresso.org, 

[22] G. Kresse and J. Purthmiiller, Gomput. Mater. Sci. 6, 15 (1996), http://www.vasp.at, 

[23] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). 

[24] V. 1. Anisimov, J. Zaanen, and O. K. Andersen, Phys. Rev. B 44, 943 (1991). 

[25] S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys, and A. P. Sutton, 
Phys. Rev. B 57, 1505 (1998), 

[26] N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12847 (1997); 1. Souza, N. Marzari, and 
D. Vanderbilt, Phys. Rev. B 65, 035109 (2001). 

[27] A. A. Mostofi, J. R. Yates, Y.-S. Lee, 1. Souza, D. Vanderbilt, and N. Marzari, Gomput. 
Phys. Gommun. 178, 685 (2008). 

[28] J. B. Goodenough, J. Solid State Chem. 3, 490 (1971). 

[29] J. R. Yates, X. Wang, D. Vanderbilt, and 1. Souza, Phys. Rev. B 75, 195121 (2007). 

[30] X. Wang, J. R. Yates, 1. Souza, and D. Vanderbilt, Phys. Rev. B 74, 195118 (2006). 


11 









[31] N. Nagaosa, J. Sinova, S. Onoda, A. H. MacDonald, and N. P. Ong, 
Rev. Mod. Phys. 82, 1539 (2010), 

[32] M. P. Lopez Sancho, J. M. Lopez Sancho, and J. Rubio, J. Phys. F 14, 1205 (1984); 
15, 851 (1985), 

[33] 1. V. Solovyev, Z. V. Pchelkina, and V. 1. Anisimov, Phys. Rev. B 75, 045110 (2007), 

[34] H. Huang, W. Duan, and Z. Liu, New J. Phys. 15, 023004 (2013), 

[35] See Supplementary Material http://link.aps.org/supplemental/xxx for more computation de¬ 
tails, the derivation of the k • p model, and a detailed analysis of the electronic structure. 

[36] T.-Y. Cai, X. Li, F. Wang, J. Sheng, J. Feng, and C.-D. Gong, arXiv: 1310.2471 (2013). 


12 





